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Abstract. We construct the spacetime in the vicinity of a general isolated, rotating, 
O ' charged black hole. The black hole is modeled as a weakly isolated horizon, and we use 

the characteristic initial value formulation of the Einstein equations with the horizon as 



£~ i ' an inner boundary. The spacetime metric and other geometric fields are expanded in a 

power series in a radial coordinate away from the horizon by solving the characteristic 
field equations in the Newman-Penrose formalism. This is the first in a series of 
papers which investigate the near horizon geometry and its physical applications using 
. the isolated horizon framework. 
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1. Introduction 



The intrinsic geometry of a classical black hole horizon in equilibrium is well understood. 
The most general treatment is provided by the framework of isolated horizons which 
^ ! allows for the possibility of a black hole with arbitrary (but time independent) intrinsic 
geometry in an otherwise dynamical spacetime [TJ EJ |3j HI |5j El [7] . Isolated horizons have 
been applied in various physical circumstances. Some illustrative examples are black 
hole thermodynamics, black hole entropy calculations in quantum gravity, numerical 
relativity, and hairy black holes. More general situations when the black grows due to in- 
falling matter and/or radiation have also been studied; see e.g. [HIE]- See [TO f [Tl j \Y2 \ [13] 
for reviews and references. 

In order to use isolated horizons in other astrophysical scenarios, it is important 
to calculate the metric in some neighborhood of the black hole. Thus it is interesting 
to solve the Einstein equations in the neighborhood of a black hole horizon. Such 
calculations have been carried out in the astropysically interesting context of a 
Schwarzschild or Kerr black hole which is tidally deformed due to its environment 
[HI [151 HS1 El HE] . Our approach here is closest to the work of Poisson and collaborators 
p~5[ H6] who study the tidal deformation of a non-rotating black hole using coordinates 
based on the past light cones originating from the horizon. These coordinates were 
originally described in jl] and are similar to the Bondi coordinates near null infinity 
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[19] . The present paper goes towards generalizing [T5J [T6] to include rotation, electric 
charge, and higher multipoles. 

We use the characteristic initial value formulation of Einstein's equations where free 
data is specified on a set of intersecting null hyper-surfaces [2Ql [211 1221 123] . Consider N 
dependent variables = 1, . . . , N) on a spacetime manifold with coodinates x a . We 
shall be concerned with hyperbolic first-order quasilinear equations of the form 



In the standard Cauchy problem, one specifies the tpj at some initial time. A solution 
is then guaranteed to be unique and to exist at least locally in time. The characteristic 
formulation considers a pair of null surfaces Ao and A/i whose intersection is a co- 
dimension-2 spacelike surface S. It turns out to be possible to specify appropriate data 
on the null surfaces and on S such that the above system of equations is well posed and 
has a unique solution, at least locally near S. 

In our case, the appropriate free data is specified on the horizon and on an outgoing 
past light cone originating from a cross section of the horizon. Such a construction in 
the context of isolated horizons was first studied by Lewandowski [23] who characterized 
the general solution of Einstein equations admitting an isolated horizon. The general 
scenario is sketched in Figure [U We consider a portion of the horizon A which is isolated, 
in the sense that no matter and/or radiation is falling into this portion of the horizon. 
For a cross-section S, the past-outgoing light cone is denoted by Af. The null generators 
of A and M are parameterized by v and r respectively; x l are coordinates on S. This 
leads to a coordinate system (v, r, x l ) which is valid till the null geodesies on M start to 
cross. The field equations are solved in a power series in r away from the horizon. This 
construction will be spelled out more precisely in the course of this paper. 

Unlike in [T5(, [16] we use the Newman-Penrose formalism [25] which, as we shall 
see, is well suited to this problem because of the central role played by null surfaces. 
Furthermore, including rotation and electric charge do not make the equations much 
more complicated in the Newman-Penrose framework. The results of this paper cannot 
yet be directly compared with [15l [16]. This would require us to start with, say 
a Kerr black hole, perturb it by considering a non-vanishing background curvature, 
and to expand the perturbation in powers of r. Perturbation theory turns out to 
be straightforward in this general set-up, at least in principle. This will be studied 
in a forthcoming paper. We make no attempt to study the issue of global existence 
of solutions, and our solutions are valid only in a neighborhood of A. There is no 
guarantee, for example, that our solutions could be extended out to an asymptotically 
flat region. Nevertheless, we expect that our construction does include most solutions 
of possible astrophysical relevance. In fact, it can be shown numerically [26] that for a 
Kerr black hole, this coordinate system extends all way out to past null infinity. It is 
thus reasonable to expect that for perturbations of Kerr, the coordinate system extends 
sufficiently far away from the horizon. 



N 




(1) 



J=l 
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Figure 1. The near horizon coordinate system. The null generators of the horizon 
are parameterized by v, surfaces of constant v are spheres. On any such sphere choose 
coordinates x 1 . The past directed outgoing light cone originating from S is J\f and r 
an afhne parameter along the null generators of A/". 

The plan for this paper is as follows. Section [2] reviews the Newman-Penrose 
formalism and summarizes the definitions and some basic properties of non-expanding 
and weakly isolated horizons. This will form the basis of the inner boundary conditions 
that will be imposed later. Section |3] sets up the near horizon coordinate system and 
gauge conditions. We start by ignoring matter fields which will be included later in 
Section [SJ Section H] summarizes the Einstein field equations and the Bianchi identities 
for vacuum spacetimes in the Newman-Penrose formalism, and Sec. |6] explicitly solves 
the field equations in powers of the radial coordinate r. Section [8] considers Einstein- 
Maxwell theory by incorporating the source free Maxwell equations. This section points 
out the specific aspects of the previous calculations which need to be modified due to 
the presence of an electromagnetic field. Finally Sec. [9] concludes with a summary and 
suggestions for further investigations. 

For simplicity, all manifolds and geometric fields shall be assumed to be smooth. 
We work in units where G — c — 1. We use the abstract index notation where lower 
case Latin letters a, b, . . . are 4-dimensional spacetime indices, and i,j, . . . denote the 
2-dimensional angular directions. We take the spacetime metric g a b to have a signature 
( — h ++). The Riemann tensor R a b c d is defined by (V a Vh — VbV a )X c = —R abd c X d 
where V a is the derivative operator compatible with g a b and X a is an arbitrary smooth 
vector field. 
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2. Basic notions 

In this section, for completeness and to set up notation, we briefly review the Newman- 
Penrose formalism and the basic definitions and properties of non-expanding and weakly 
isolated horizons. 

2.1. The Newman- Penrose formalism 

The Newman- Penrose formalism [25l [27] is a tetrad formalism where the tetrad elements 
are null vectors, which makes it especially well suited for studying null surfaces. See 
[28| |29| 123] for pedagogical treatments (note that these references take the spacetime 

metric to have a signature of (H ) which is different from ours). Start with a null 

tetrad (£, n, m, m) wher £ and n are real null vectors, and m is a complex null vector 
and m its complex conjugate. The tetrad is such that £ ■ n = — 1, m ■ m = 1, with all 
other inner products vanishing. The spacetime metric is thus given by 

9ab = ~Ln b ~ n a £ b + m a m b + m a m b . (2) 

Directional derivatives along the basis vectors are denoted as 

D:=£ a V a , A:=n a V a , 5:=m a V a , 5:=m a W a . (3) 

(The symbol A is used for both the directional derivative along n a and for the isolated 
horizon). The components of the connection are encoded in 12 complex scalars, the spin 
coefficients, defined via the directional derivatives of the tetrad vectors: 



D£ 


= (e + e)£ — Rm — nm , 


(4a) 


Dn 


= — (e + e)n + rem + fern , 


(4b) 


Dm 


= ni — Kn + (e — e)m , 


(4c) 


A£ 


= (7 + j)£ — fm — rm , 


(4d) 


An 


= —(7 + j)n + vm + vm , 


(4e) 


Am 


= v£ — rn + (7 — 7)m , 


(4f) 


5£ 


= (a + (3)£ — pm — am , 


(4g) 


5n 


= —(a + j3)n + fim + Am , 


(4h) 


Sm 


= \£ — an + (/? — a)m , 


(4i) 


5m 


= fi£ — pn + (a — (3)m . 


(4j) 



V 



Many of the spin coefficients have a transparent geometric interpretation. Some 
important ones for us are: the real parts of p and p are the expansion of £ and n 
respectively; a and A are the shears of £ and n respectively; the vanishing of k and 
v implies that £ and n are respectively geodesic; e + e and 7 + 7 are respectively the 
accelerations of £ and n, a — (3 yields the connection in the m-m plane and thus the 
curvature of the manifold spanned by m-m. 
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Since the null tetrad is typically not a coordinate basis, the above definitions of the 
spin coefficients lead to non-trivial commutation relations: 



(AD - DA)f 
(5D - D6)f 
\ (6A-A6)f 
(66 - 66) f 



(e + e) Af + ( 7 + j)Df - (f + ir)6f - (r + rr)6f , (5a) 

(a + (3- n)Df + nAf - (p + e - e)6f - a5f , (5b) 

-uDf + (r-a- (3)Af + (jj, - 7 + 7)6 f + XSf , (5c) 

ifi - y)Df + (p- p)Af + (a- 13)6 f -(a- (3)6 f . (5d) 



The Weyl tensor C a bcd breaks down into 5 complex scalars 



^0 = C abcd £ a m b £ c m d , ^ = C abcd £ a m b £ c n d , ^ 2 = C abcd £ a m b m c n d . 
^3 = C abcd tn h m c n d , # 4 = C abcd m a n b m c n d . 



(6a) 
(6b) 



Similarly, the Ricci tensor is decomposed into 4 real and 3 complex scalars <J>jj: 

= l -R ab £ a £ b , $ n = \R ab (tn b + m a m b ) , $ 22 = ^R ab n a n b , A = ^ 

1 1 1 

= -R ab £ a m b , $02 = -Rabm a m b , $12 = -R ab m a n b , = $, 




(7a) 
(7b) 



The 6 components of the Maxwell 2-form F ab are written in terms of 3 complex scalars 



-F ab £ a m b , 



^F ab (n a £ b + m a m b 



F ab n a m 



(8) 



The source free Maxwell equations dF = and d*F = are written as 4 complex scalar 
equations 



Acj) 
A0i 



6 = (7r — 2a)(f) + 2p<p l - K(f) 2 , 
b x = -X(p + 27T0! + (p - 2e)(p 2 

f»l = (27 - /i)0 o - 2T01 + (T02 , 



z/0o - 2/i0! + (2(3 - r)0 2 . 



(9a) 
(9b) 
(9c) 
(9d) 



The stress energy tensor for the Maxwell field is given by 



T ab — — [ F ac F b 



-^9abF cd F 



cd 



(10) 



which is seen to be trace-free. Using the Einstein equations, we see that for Einstein- 
Maxwell theory, the Ricci tensor components defined in Eqs. ([7]) have a simple expression 
in terms of the cbf. 



2^ (i,j = 0,1,2), A = 



The relation between the spin coefficients and the curvature components lead to the 
so called Newman- Penrose field equations which are a set of 16 complex first order 
differential equations. The Bianchi identities are written explicitly as 8 complex and 3 
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real equations. See J2S] 122] 123] for the full set of field equations and Bianchi identities. 
We shall later write these equations after imposing gauge and coordinate conditions. 

It will also be useful to use the notion of spin weights and the 9 operator for 
derivatives in the m-rh plane (which will be angular derivatives in our case). A tensor 
X projected on the m-rh plane is said to have spin weight s if under a spin rotation 
m — > e^m, it transforms as X — > e ts ^X. Thus, m a itself has spin weight +1 while m a 
has weight —1. For a scalar X = m ai ■ ■ ■ m ap rh bl ■ ■ ■ rh bq X ai ...\ }q , i.e. it has p contractions 
with m and q with fh, then X has spin weight s = p — q. For example, the Weyl tensor 
component has spin weight 2 — k. Similarly, the Maxwell field component <pk has 
weight 1 — k. 

The S and 5 operators are defined as 

clX = m ai ■ ■ ■ m a "rh bl ■ ■ ■ rh b «5X ai ... bq , (12) 
QX = m ai ■ ■ ■ m a ^m bl ■ ■ ■ m b «5X ai ... bq . (13) 

From Eqs. OH]) and (4jJ, after projecting on to the m-rh plane, we get 

Sm a = {(3 - a)m a , 5m a = {a - P)m a . (14) 

A short calculation shows that 

QX = 5X + s(a - 0)X , dX = 5X- s(a - j3)X . (15) 

It is clear that 5 and 5 act as spin raising and lowering operators. See [30J for further 
properties of the 9 operator and its connection to representations of the rotation group. 
The transformations of the null tetrad which preserve the metric are 

(i) the boosts: 

/ — > Al , n — ?• A~ 1 n , m — > m , (16) 

(ii) the spin transformations in the m — m plane: 

m-J-e^m, £—>■£, n n , (17) 

(iii) the null rotations around £: 

m— >m + a£, n — >■ n + am + am + \a\ 2 £ , (18) 

(iv) the null rotations around n (obtained by interchanging £ and n in Eq. (fl8]) ). 
Again, we refer to [28] 129] [23] for a more complete discussion. 



2.2. Non-expanding and weakly isolated horizons 

A black hole in equilibrium with its surroundings is modeled quasi-locally as an isolated 
horizon. The basic geometrical object is a smooth 3-dimensional null surface A (which 
shall be the black hole horizon) in a Lorentzian spacetime (Ai, g a b)- If £ is any null 
normal of A, then it must be geodesic so that £ a V ' a £ b = H^£ b ; the acceleration Km is 
the surface gravity associated with £ a . We shall consider only non-extremal horizons 
here, i.e. we shall always have non-vanishing k. The spacetime metric g a b induces a 
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degenerate metric on A which we denote g a & which has signature (0 + +). Its (non- 
unique) inverse will be denoted q ab . There is a volume element 2 e a t > on A, satisfying 
£ a2 £ a b = 0, which measures the area of spacelike cross-sections of A. 

A smooth 3-dimensional null surface A is said to be a non- expanding horizon if: 

• A has topology S 2 x R. 

• The expansion 0«) := q ab, V a £b of any null normal £ a of A vanishes. 

• The Einstein field equations hold at A, and the matter stress-energy tensor T a j, is 
such that for any future directed null-normal £ a , —T££ b is future causal. 

We shall consider only null tetrads adapted to A such that, at the horizon, £ a coincides 
with a null- normal to A. We shall also consider a foliation of the horizon by spacelike 
spheres S v with v a coordinate on the horizon which is also an affine parameter along 
i: Civ = 1. Null rotations about £ a correspond to changing the foliation. 

This deceptively simple definition of a non-expanding horizon leads to a number 
of important results which we state here without proof (though some of these will be 
rederived later): 

(i) The Weyl tensor components \l/o and vanish on the horizon. This implies that 
^2 is an invariant on A as long as the null-tetrad is adapted to the horizon; it is 
automatically invariant under boosts and spin rotations (it has spin weight 0), and 
it is invariant under null rotations around £ because \l/o an d ^i vanish. Similarly, 
the Maxwell field component 0o vanishes on the horizon, and 0i is invariant on A. 
Both \&2 and 4>i are also time independent on the horizon. 

(ii) For a general null sub-manifold, there is no unique derivative operator compatible 
with the metric, and the pull back of the spacetime derivative operator V a does not 
necessarily induce a connection on the hypersurface. For non-expanding horizons 
however, the spacetime connection does induce a unique derivative operator 
compatible with q^. Furthermore, there exists a 1-form u a such that, for any 
vector field X a tangent to A, 



The 1-form oo a plays a fundamental role in what follows. The pullback of oj a to the 
cross-sections S will be denoted u a . 

(iii) The surface gravity of £ is 



The curl and divergence of oo carry important physical information. The curl is 
related to the imaginary part of the Weyl tensor on the horizon 



X a V a t = X a uJ b . 



(19) 



a 



(20) 



du = Im [# 2 ] 2 e . 
and its divergence specifies the foliation of A by spheres [6]. 



(21) 
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(iv) The horizon angular momentum is well defined in the case when there is an axial 
symmetry (p a on A which preserves q a b, u a and the electromagnetic field on the 
horizon [5]. The angular momentum is given by 

J = -^£ / Im ^ 2 ] h + ^f s ^ Im [0i] 2 e , (22) 

where / and g are respectively defined via tp a2 e a b = dbf and (p a * F a b = dbg. 
Similarly, the electric and magnetic charges of the horizon are defined respectively 

as 

Q = — I Re 2 e , P = — t Im [fa] 2 e . (23) 



27T Js 27T j 5 

Hamiltonian methods provide a suitable notion of horizon mass [5] : 

M=^vW + W+^, (24) 
where R is the area radius of the horizon so that if A is the area of the horizon cross- 



sections then R := ^A/An. In the case when the horizon is not exactly symmetric, 
one could attempt to find an approximate symmetry to replace (p a in the above 
equations [311 [32j [33]. It is also possible to define source multipole moments for an 
axisymmetric charged isolated horizon [34]. 

We need to strengthen the conditions of a non-expanding horizon for various physical 
situations. The minimum extra condition required for black hole thermodynamics and 
to have a well defined action principle with A as an inner boundary of a portion of 
spacetime, is formulated as a weakly isolated horizon [3]. This is to choose an equivalence 
class of null normals [£}, each related to the other by a constant re-scaling £' = ci, such 
that 

C#j a = . (25) 

This can be shown to be equivalent to the zeroth law, i.e. Km = £ a u a is constant on the 
horizon. Note that under a re-scaling £ a — > f£ a , u a transforms as co a — > u a + D a \nf so 
that it is invariant under constant rescalings. This condition is sufficient to ensure a well 
defined action principle and to lead to a sensible notion of horizon mass and spin. Any 
non-expanding horizon can be made into a weakly isolated horizon by suitably scaling 
the null generators. Thus, the restriction to weakly isolated horizons is not a genuine 
physical restriction. 



3. The near-horizon coordinate system and null-tetrad 

Let us now assume that the vacuum Einstein equations hold in a neighborhood of the 
horizon A. We will consider electromagnetic fields later in Sec. [H Following [1] we 
introduce a coordinate system and null tetrad in the vicinity of A analogous to the 
Bondi coordinates near null infinity. See Fig. [U Choose a particular null normal £ a on 
A. Let v be the affine parameter along £ a so that £ a V ' a v = 1. Let S v denote the spheres 
of constant v. Introduce coordinates x % (i = 2, 3) on any one S v (call this sphere So) 
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and require them to be constant along £ a : fV tt i ! = 0; this leads to a coordinate system 
(v,x l ) on A. Let n a be a future directed inward pointing null vector orthogonal to the 
S v and normalized such that I ■ n = — 1. Extend n a off A geodesically, with r being an 
affine parameter along — n a ; set r = at A. This yields a family of null surfaces M v 
parameterized by v and orthogonal to the spheres S v . Set (v,x l ) to be constant along 
the integral curves of n a to obtain a coordinate system (v, r, x l ) in a neighborhood of 
A. Choose a complex null vector m a tangent to S . Lie drag m a along £ a : 

C e m a = on A. (26) 

We thus obtain a null tetrad (£, n, m, m) on A. Finally, parallel transport i and m along 
— n a to obtain a null tetrad in the neighborhood of A. This construction is fixed up to 
the choice of the x l and m a on an initial cross-section Sq. We are allowed to perform 
an arbitrary spin transformation m — > e^m on S . 

It is easy to write the elements null tetrad in these (v, r, x l ) coordinates. Let us 
start with n a and n a . We have the family of null surfaces M v parameterized by v; n a 
is normal to the J\f v , and r is an affine parameter along — n a . This implies that we can 
choose 

n a = -d a v and n a W a := A = ~ . (27) 

Or 

To satisfy the inner-product relations, the other basis vectors must be of the form: 

rV o: = D = | + 4 + A-<J-, m «V^ = 4 + fJ-. (28) 

The frame functions U, X 1 are real while f2, ^ are complex. We want £ a to be a null 
normal of A so that the null tetrad is adapted to the horizon. Since d v is tangent to 
the null generators of A, this clearly requires that U,X % must vanish on the horizon. 
Similarly, we want m a to be tangent to the spheres S v at the horizon, so Q should also 
vanish on A. Thus, U, X l ,Q are all 0(r) functions. 

We expand the spin coefficients, Weyl tensor components and the directional 
derivatives in a power series in r away from the horizon 

X = X^+rX^ + -r 2 XW... . (29) 

Thus, for example we will have 

* k = *f + r*™ + + .... (30) 

The same notation will be used for the frame fields U and f2. However, for the frame 
fields X 1 and £ l , we will write 

e = e { o)+re (l) + \r% ) + ... . (31) 

To avoid clutter we shall not be completely consistent with this notation. Thus, we shall 
not use any index for the spin weighted angular derivative 9; it is to be understood that 
9 always refers to in this paper. Similarly, where we don't expect any confusion, we 
shall often drop the index on the directional derivatives such as D and 5; in this case, 
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unless mentioned otherwise, the relevant order of the operator is the same as the order 
of the operand. For example, S^f^ refers to . 

The variables we need to solve for are the frame fields U, X l ,Q,l; 1 , the 12 spin 
coefficients and the Weyl tensor components \l/fc. The equations are the commutations 
relations Eqs. (128fl . the 16 field equations and 8 of the 11 Bianchi identities. Three of 
the Bianchi identities involve only the Ricci tensor, thus they will need to be considered 
when matter fields are present. In Sec. [8] we will consider the Maxwell equations and the 
3 additional Bianchi identities as well. All these are first order differential equations and 
each of these sets of equations has possibly three kinds of equations: evolution equations 
which involve derivatives along v , i.e. D, and do not contain any radial derivatives A, 
equations which contain only purely angular derivatives 5 and 5, and finally the radial 
equations involving A. In order to integrate these equations, we proceed as follows. 
We start with suitable data on some initial cross section So of the horizon and use 
the non-radial equations to propagate them at all points of the horizon. Starting with 
this horizon data and the appropriate data on the past light cone Ao containing So, 
the radial equations then yield the first radial derivatives and, iteratively, all successive 
higher derivatives as well. At each step, we will need to expand the non-radial equations 
in powers of r, and ensure that we have consistency order-by-order in r. We now carry 
out this procedure in detail and spell out the free data. 

We start with the conditions on the spin coefficients. Since n a is an affinely 
parameterized geodesic, and I and m are parallel propagated along n a , we have 
An = A£ = Am = 0. From Eqs. did]) . (He)) and flifj), this leads to 

<y = T = V = 0. (32) 

Imposing Eq. (|32|) in the commutation relations ()5]) leads to 

(A-D — DA)f =( e + e)Af-ir5f-7r5f, (33a) 

(SD-DS)f = (a + (3 -it)Df + kAJ - (p + e-e)8f - aSf , (33b) 

(5A-A8)f — —(a + /3)Af + p5f + X8f , (33c) 

(55 - 55)f =(ji-ii)Df + (j}-p)Af+(a-P)6f-(a-P)6f. (33d) 

Substituting f = v yields 

ii = a + (3 , /i = /2. (34) 

Eqs. fl32|) and (134"]) are the basic conditions on the spin coefficients which hold at 
all spacetime points where the coordinate system is valid. In addition, we have the 
boundary condition at the horizon, namely that £ is geodesic, normal to a smooth 
surface, and expansion free. This implies 

p(o) = /c (o) = . (35) 

The subscript ^ indicates that these are the values at the horizon, i.e. at r = 0. From 
the Raychaudhuri equation applied to the null generators of the horizon, we get 

Dp® - 6k® = (p (0) ) 2 + + (e<°> + eW)p (0) - 2a (( V 0) (36) 
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(this is also one of the Newman- Penrose field equations considered later). This leads to 
cr(°) = 0. In addition, we can make a spin transformation on Sq and set 



From Eq. (I33bl) and the above conditions on the spin coefficients at the horizon, this 
also ensures that Lim a = at the horizon. Finally, we assume that we have scaled the 
null normal appropriately (i.e. we have chosen the coordinate v and the foliation S v ) so 
that the horizon is weakly isolated: L(UJ a = 0. In terms of the Newman- Penrose spin 
coefficients, u a is written as 



Among the free functions that we need to specify, we can choose a constant for the 
surface gravity (perhaps determined by the mass and spin (5j[3]), and a function tt^ on 
the initial cross-section So. 

4. The reduced system of equations 

With the conditions Eqs. fl32|) and (1341) on the spin coefficients at hand, we now impose 
them in the commutator relations, the field equations and the Bianchi identities. The 
functions U, X\ Q, £ l are determined by the commutation relations ( |33l) by substituting, 
in turn, r and x l for /, and imposing Eqs. (|32|) and (I34|) on the spin coefficients. First 
the radial derivatives: 



e (0) _ e -(o) = _ 



(37) 



w a = ~n b W a t = (e + e)n a + 7r (0) m a + 7r (0) m, 
Thus, A will be a weakly isolated horizon if we choose 
k ■= e (°) + g(°) = constant , Dn {0) = . 



a 



(38) 



(39) 



< 



AC 



AU 
AX i 

An 




(40a) 
(40b) 
(40c) 
(40d) 



v 



The propagation equations along v are 




(41a) 
(41b) 



Let us now turn to the field equations. After imposing Eqs. (132]) and ( |34|) on the 
spin coefficients and ignoring matter terms for now, the field equations involving radial 
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derivatives are: 



AA = -2X/J, - * 4 , 
Ap = -p 2 -\X\ 2 , 
Ap = -fip - crA - * 2 , 
Aa = —pa — Xp . 
Ak, = —Tip — 7ro" — * x , 
Ae = -Txct - ir(3 - * 2 , 
An = —TTfl — ft A — *3 , 

A/3 = — p(3 — aX , 

Aa = -/3A - pa - * 3 . 



The time evolution equations become: 



Dp — 5n = p 2 + \a\ 2 + (e + e)p — 2an , 

Da-5K=(p + p + 2e)a - 2(3k + * , 

Da — 5e = (p — e)a + (3a — (3e — kX + (e + p)n , 

D(3 — 5e = (a + n)a + (p — e)/3 — pn — (a — n)e + *] 

DX-H = (p- 2e)A + o-/i + 2avr , 

Dp - 5tt = (p - 2e)p + o-A + 2/3tt + * 2 . 



The angular field equations are: 



5p — 5a = np — (3a — (3)a — *i , 

5a - 5(3 = pp - \a + \a\ 2 + \(3\ 2 - 2a(3 - * 2 , 

5\ - 5p = Tip + (a - 3/3)A - * 3 . 



(42a) 
(42b) 
(42c) 
(42d) 
(42e) 
(42f) 
(42g) 
(42h) 
(42i) 



(43a) 
(43b) 
(43c) 
(43d) 
(43e) 
(43f) 



(44a) 
(44b) 
(44c) 



Finally, we have the Bianchi identities which, in the NP formalism, are written as a 
set of nine complex and two real equations; in the absence of matter, only 8 complex 
equations survive. The radial Bianchi identities reduce to: 



A* - #*i 
A* x - 5* 2 

< A* 2 - 5* 3 
A* 3 - <5* 4 



-/i*o - 2/3*1 + 3fr* 2 , 

-2pm l + 2(7*3 , 
-3/i* 2 + 2/3*3 + <^*4 , 

-4/i* 3 + 4/3*4 • 



(45a) 
(45b) 
(45c) 
(45d) 



Note that there is no equation for the radial derivative of *4. Among all the fields that 
we are solving for, this is in fact the only one for which this happens. This means that 
*4 (in this case, its radial derivatives) is the free data that must be specified on the null 
cone A/q originating from S . 
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Finally, we have the components of the Bianchi equations for evolution of the Weyl 
tensor components: 





-6V = 


(tt - 4a)^ + 2(2p + e)*j - 3/c* 2 , 


(46a) 




= 


-A^ - 2(7r - a)*i + 3p^ 2 - 2k* 3 , 


(46b) 




- 5^2 = 


-2A*i + 3n^ 2 + 2(p - e)^ 3 - , 


(46c) 




-5^ 3 = 


-3A^ 2 + 2(a + 2vr)^ 3 + (p - 4e)^ 4 • 


(46d) 



V 



Having written down the full set of differential equations, we are now ready to solve 
them locally near the horizon order-by-order in r. 



5. The intrinsic horizon geometry 

We first consider the intrinsic geometry of the horizon contained in the time evolution 
equations at the horizon, i.e. the O(r ) part of Eqs. ( 1431) . This has been studied 
elsewhere in great detail [6j [3j [7] , so we shall be brief. We already have that £ a is geodesic 
so that k® = 0, and the expansion and twist of i vanish, which implies p® = 0. The 
Eqs. ( 14"3|) in turn lead to 

' ^(0) = o , ^ 0) = , (47a) 

Da® - 6e® = e®a® - P®e® + e®ir® , (47b) 
D p(o) _ Se (o) = (q (o) + ^(0)^(0) + e (o )/3 (o) _ ( -(o) _ - ( o ))e (o) + ^(0) ? (47c) 

DX® - 8tt® = -2e^A^ + 2a(°)7r(°) , (47d) 
Dp® - Sir® = -2e®p® + 2{3®tt® + *f . (47e) 



Adding Eq. (I47bj) to the complex conjugate of Eq. (|47c| . and using the conditions on 



the spin coefficients leads to 

Da® = , £>/3 (0) = , = . (48) 

The remaining two equations tell us about the time dependence of the expansion and 
shear of n at the horizon: 

DX® + RX ® = ^(o) + 2a (o) 7r (o) j (49) 

ZV°) + «p(°) = 5tt(°) + 2/3<V°> + *f • (50) 

Using the 9 operator for the angular derivatives, and noting that 7r has spin weight — 1, 
this is written more cleanly as 

DX® + RX® = 3ir® + (vr^) 2 , (51a) 

Dp<°> + £p (0) = cV°) + \tt®\ 2 + ttf . (51b) 
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These two equations can be shown to be identical to Eq. (3.9) of [6] with the matter 
terms set to zero. 

Now consider the angular field equations (jUJ) at the horizon. The first of these 
leads to \E f ^°- ) = which we already knew. Eq. fl44c[) gives 

#(°) = 8^(0) + ^(0)^(0) _ g A (0) _ 7r (0) A (0) (52) 

The real and imaginary parts of Eq. (144b j) give respectively 

-2Re^ 0) = 5a + 5a - 2\a\ 2 , a := a (0) - ^ (0) , (53) 

-2zlm*f = 5vr (0) - e7r (0) . (54) 

Since a = — determines the connection on the horizon cross-section, it's 
derivative must be connected with the scalar 2-curvature of the cross-section. In fact, 
it is not hard to show that the first equation is equivalent to 

2 TZ = -4Re^ 0) (55) 

The second equation relates the curl of Co with the imaginary part of ^2 and thus 
reproduces Eq. ( 12~T|) . 

Now the Bianchi identities Eqs. fj46|) lead to 



D^f = , (56a) 
+ (°> = tf* (°> + 3n^¥ 2 0) , (56b) 
Dttf + 2K*f = N 0) + 57r(°)*f - 3A«%< 0) . (56c) 



The second of these can be shown to follow from Eq. ( 1521) . so it doesn't yield any new 
information, but the last equation gives the time evolution of ^4 at the horizon. 

To summarize, here is the free data required at the horizon. We start with a 
choice of null generator t a and a constant surface gravity k := e*- -* + e*- -*, and an affine 
parameter v. Surfaces of constant v are spheres S v , and we choose a basis m, m tangent 
to the S v and Lie dragged along £ a . Then, on any one of the spheres, say So, choose 
the transversal expansion and shear \(°\ the spin coefficient the connection on 
So'- <y^ — and the transverse gravitational radiation ^4 . The field equations then 
show that 7P ', cr°) — /3 < -°- ) are time independent. The transversal shear and expansion 
\W(v),fj,V>{v) satisfy Eqs. (EI]) which means that their time evolution is: 



^)( v ) = ^°)(0) e -' m + - [cV°) + (vrW) 2 ! (1 - e- kv ) , (57a) 

K 



\W(v) = \ (0) (0)e~ kv + I [Stt^ + Itt^I 2 + ¥ 2 0) 1 (1 - e~~ KV ) . (57b) 

The Weyl tensor components ^0 ,^1 vanish. \E f 2°' > is time independent, it's real 
part is determined by — according to Eq. ( ISTj) . and its imaginary part by 
according to Eq. ( 154]) ; \&g is determined from Eq. (152]) and its time 



evolution is of the same form as for A*- -*, From Eq. (146bj) we get 



yg>(v) = #f (0) e - Rv + 4 \m { 2 0) + 37r (0) ^f 1 (1 - e~' KV ) . (58) 



K 1 
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Finally, can be freely specified on Sq and it evolves in time according to Eq. (I56cl) ; 
this leads to a time dependence of the form 



*i 0) (0) + ^(l-e 



A 



-2kv 



B 
+ — e 

K 



where the time independent angular functions A and B are defined via 
(8 + 5tt(°)) ttf - 3A^ 0) = A + Be~™ , 



A 



1 



(8 + 57r (0) )(B + 37r (0) ; 



B = (8 + 57r(°))v|4 0) (0) - 3A (0) v|/f - A . 



3(9 

(0) 



ITT 



(0) 



(0)1 



(0) 



(0) 



(59) 

(60) 
(61) 
(62) 



We note that the time evolution on an extremal horizon (k = 0) would be very different; 
see e.g. [6j [7]. As we shall see in the next section, not only can we specify \l/ 4 freely on 
Sq, but also all of its radial derivatives. The free data on the transversal null surface 
Aq is in fact \&4. 



6. The radial equations 

With the intrinsic horizon geometry understood, we turn our attention to the radial 
equations (H2l) and (H5l) which determine all radial derivatives for all the non-zero spin 
coefficients and Weyl tensor components. The first radial derivatives are easily obtained 
by substituting the horizon values on the right hand sides of Eqs. ( 1421) and ( |45l) : 



(i) 



(i) 
(o) 







P 

e (D + e -(D 



2|7r (0) r + 2Re[# 



(0)1 



,(i) 



— e 



(i) 



(0) 



A (1) = + * 4 

/U (D = ( /U (o))2 + |a(°)| 2 



7f(o)( a (o) - ^(°)) - tt^c^ - + ¥^ - * 



(0) 



7T 



a: 



(i) . 
(i) 



(i) 



-0, 

-- -m 



+ vr (o) A (o) + vl> 
^(o)) + A (o) (/3 (o)_ 



(0) 
3 J 

- # 



(0) 



(5 + 7r(°))^ 0) + 3/i( ^ 
(5 + 2vr(°))^i 0) + V°)^ 



(o) 



(0) 
3 



(63a) 
(63b) 
(63c) 
(63d) 
(63e) 
(63f) 
(63g) 
(63h) 
(63i) 
(63j) 
(63k) 
(631) 



We can continue this process by applying the radial derivative A to Eqs. (jl2"j) once again, 
and substituting the first derivatives from above on the right hand side. This leads to 
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expressions for the second derivatives in terms of the horizon data: 



C0 = -(3-fr<°))4f 



(2) = \ (0)^(0) 



p (2) 

A (2) 



(0) 



*(0)W°) 



(0) 



6( M W) 2 AW +2A^ |A^| + 2//<% 
2(^)3 + 6/i (o) | A (0)| 2 + A (o)^(o) + A (o)^ 



(1) 



r (0) 



(0) 
3 



7r( 2 > = 27r(°)((/i( )) 2 + |At°>P + 4vr(°)A( V( ) + 5^°>* 
+A (o)^(o) _ (a + ^ (o) )xI/ (o) ^ 

a (2) _ ^(2) = 2(a ( 0) _ / g(o))((^(o))2 + | A (o)| 2 ) + 4/ /o) A (o) (/3 (o) - «(")) 



0* 



(0) _ fi (0) 



) + 5^°>ttf - A^f - (3 + 2tt(°))^ 



(64a) 
(64b) 
(64c) 
(64d) 
(64e) 

(64f) 
(64g) 



Finally we can also investigate the non-radial equations to 0(r). For this, we need the 
directional derivatives to 0(r) which we have already calculated: 



A - - — 

Or ' 

D - I + ^ + ^Wr + + + ° (r2) 



><®dx 



(o) dx i 



dr 



r 



(65) 
(66) 
(67) 



We substitute these on the left hand sides of Eqs. ( HIT) . ( 1431) and ( 1441) . and also expand 
the right hand sides in powers of r. Equating powers of r on both sides yields differential 
equations at each order. The O(r ) terms have already been studied in Sec.[5]since these 
deal with the intrinsic horizon geometry. The 0(r) equations are straightforward in some 
cases. For example, Eq. (I43bj) is trivial because k, a, \l/ are ah 0(r 2 ) quantities as we 
have already seen in Eqs. f lBHj) . Similarly, Eq. (143 ap is also easy. Using p = rmf + 0{r 2 ) 
and the vanishing of k, o to 0(r), we get 



— * 

dv 



(0) 







(6£ 



which we already knew. The remaining equations are similarly straightforward, but a 
lot more tedious. In the end, it turns out that Eqs. (|4T|) . ( [43]) and ( 14"4")) do not yield any 
new information at 0{r). 



7. The expansion of the metric 

To expand the metric in powers of r, we start with the frame fields defined in Eqs. (p7|) 
and (I28p . and the radial frame equations (140!) derived from the commutation relations. 
The strategy is the same as for the spin coefficients. Eqs. (1401) give us the first radial 
derivatives by substituting the horizon values on the right hand side, and taking higher 
derivatives leads to the higher order terms. The calculations are straightforward and 
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U = rk + r 2 {2 |tt (0) f + Re[*f ]) + C(r 3 ) , 
Q = rTfW + r 2 ^(°)7f(°) + A< M°> + + G(r 3 ) , 

= O^ 0) + ng 0) + 0(r 3 ) , 

e= [l + r/i (0) +r 2 ((^)) 2 + |A(°)| 2 )]^o) 
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(69a) 

(69b) 
(69c) 



+ 



(0) 



V ' 2 

The contravariant metric is seen to be given in terms of the frame fields as follows: 



(69d) 



g rr = 2(U + \n\ 2 ) , 
9 vr = l, 

g ri = x ,i + nc + nc 

g^ = ee+ee- 



An explicit calculation yields 



9' 

9 l 



2«r + 2r 2 (3|tt (0) | 2 + Re[*f ]) + £>(r 3 ) . 
2rvr (o) + r 2 f 2^(0)^(0) + 2A (o )7r (o) + 1^(0) A + 0{? 



(70a) 
(70b) 
(70c) 
(70d) 



(71) 



f(0) 



+ 



2r7r (o) + r 2 / 2^(0)^(0) + 2^(0)^(0) + l^(o)A + 0(r . 3) 



(72) 



;(o) 



g* = (1 + 2r^°) + 3r 2 ((/i(°)) 2 + |A(°)| 2 )) + ^ 0) ) 



+ 
+ 



rA^+r 2 ( 3/i(°)A(°) + ^ 



(0) 



rA(°)+r 2 (3/i( )A(°) + ^ 



(o) 

4 



2 ^(o)C(o) 
2^o)§o) + 0(i 



(73) 



The null co-tetrad, i.e. the dual basis for the null tetrad found above can also be 
calculated easily up to 0(r 2 ) : 

n = -dv , (74) 
i = dr - (kr + Re[^ 0) ]r 2 ) dv - (V V + ^ 0) r 2 ) tfW - + ^3°V 2 ) g°W , (75) 

m = - (V°V + r 2 ^j dv + (1 - //°V)f <°W - (Vv + ^fr 2 ^) g°W . (76) 

Here £^ are defined by the relations £i°^(o) = and ^-^l^ = 1; it will be convenient 
to set mi°^ := ^d a x l . In powers of r, the metric is: 

g ab = -2l [a n b) + 2m {a m b) = g^ + + ]^r 2 + ■■■ , (77) 
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where 

g^ = 2d (a rd b) v + 2m^fn^ ) , (78) 

g$ = - (2kd {a vd b) v + A7T^m^d b) v + An^fhfjd^v 

+ A^m^m^ + 2X^m^m^ + 2X^m^m^) , ^ 

g^=A{\^-Re[^])d {a vd b) v 

+ 4 (>V ) + A<%< > - mg^i; 

+ 4 f// (0) 7f (0) + A<°M°> - *<°>) mgffyi; (80) 
+ 4 ((m (0) ) 2 + |A (0) | 2 ) mjM? 

+ (4//°)A(°) - 2*f) ro£>rng> + (V°>A<M - 2*f) mj^mg 5 . 

We could easily obtain the third order metric by using Eqs. ( |64|) . but we shall not do so 
here. 



8. Inclusion of electromagnetic fields 

Let us now go beyond vacuum spacetimes and include electromagnetic fields. As usual, 
we expand the Maxwell field components <pk in powers of r: 

A = <f>k ] + r <f>k ] + + ■■■ > where k = 0,1,2. (81) 

Consider first the intrinsic horizon geometry. We start with the behavior of the Maxwell 
field components at the horizon. The Raychaudhuri equation for the null generators, 
after imposing the conditions of Eqs. ( )32l) and ( J34l is just one of the Newman- Penrose 
field equation (Eq. l43aH that we have used earlier. In the presence of a Maxwell field, 
that equation becomes 

Dp-5n = p 2 + H 2 + (e + e)p - 2an + 2\<p \ 2 . (82) 

On the horizon, by definition p = k = 0. Thus from the above equation we deduce that 
we still have a = as before, and in addition we also get 



. (83) 



This implies that 

*® = *® = *® = *® = *® = 0. (84) 

Furthermore, we also have A = everywhere since the trace of T ab vanishes. Imposing 
these conditions to determine the time evolution of the geometrical fields on the horizon, 
the two non-radial Maxwell equations ( l9aj) and ( l9bl) become 



D0i o) = , (85a) 
D$ + = M 0) + 2vr(°)0S o) . (85b) 
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Thus, <fii is time independent, the charges defined in Eq. ( l23|) are constant. Since 0^ 
vanishes, these charges are in fact independent of which cross-section of A is used in 
the integral. Just like ^3 and ^4 at the horizon, 2 has a time dependence. 

With these simplifications at hand, it turns out that all of Eqs. (I47j) remain 
unchanged. The only change in the intrinsic geometry occurs in two of the angular 
equations ( )44l) . The first change is an extra contribution in Eq. ( 1531) : 

-2Re[^ 0) ]+4|^ 0) | 2 = 5a + 5a-2\a\ 2 -4Re[^ 0) ] + 8\(j)f ) \ 2 = 2 ft.(86) 

The second effect is in Eq. (1521) which determines 

«rW = g M W + 7r(V 0) - 3A W - 7f(°)A(°) + 20f0j o) . (87) 

The remaining Maxwell equations ( J9cl) and ( I9dl) determine the radial derivatives of 
0o and 0i. As usual, we impose Eqs. (1321) and ( 1341) to obtain the radial equations for O 
and 0i : 

A0 O = ^01 - W>o + c02 , (88a) 
A0i = 502 - 2/^0! + 2/302 • (88b) 



Note that there is no equation for the radial derivative of 02- Thus, just like ^4, we 
need to specify -02 on the past light cone Ao as part of the free data. As for the other 
radial equations, we can iterate this to obtain an expansion for O and 0i. To first order 
we need to simply evaluate the right hand side of these equations at the horizon: 

/>« = -tyf) , (89a) 
6 W = + 2/U (°)0j°) + 2/?(°)0f = (-5 + 7r(°))0( o) + 2/i(°)0S o) . (89b) 

Let us turn now to the Bianchi identities. In the vacuum case, we had the Equations 
( I4l)|) and (146 p . These equations give here as well the radial and time derivatives of the 
Weyl tensor components with additional terms for the matter contributions. For the 
time derivatives at the horizon, we only need to worry about (since #J 0) has been 
determined above and the others are time independent): 

+ 2«*?> = (3 + 5nW)¥ 3 0) + 20f )) (B + 

_ 3A (O) vI/ (O)_ 4A (O)| (O) | 2 > W 

The radial equations for the Weyl tensor components from the Bianchi identities are not 
very enlightening, and we shall just write the first derivatives evaluated at the horizon: 

*f ) = 0, (91a) 

*W = _ 20f ) 50i o) , (91b) 

= _(g + «(o))y£) _ 20S O) (5 + vr(°))0( o) + 20( o) 50i o) + 3/xtf< 0) + 4/i|0f ) | 2 , (91c) 
vD« = _(a + 27r<°>)* < 0) + 4^f + + 2<g 0) (3 - 2^)^ 

-4A(°)0f ) 0f +47r(°)|0f| 2 (91d) 
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In the presence of matter fields, there are three (two real and one complex) additional 
Bianchi identities containing only Ricci tensor components. Let us consider these 
equations at the horizon after imposing our conditions on the spin coefficients and the 
vanishing of 0o- The first of these equations yields simply that 0q ^ is time independent 
(which we already knew): 



(0) 

11 



=> Dcf) 



(0) 

i 



0. 



The second determines the time evolution of $ 



(0) 
12 



(o) 

12 



The third yields the time dependence of $22: 



(0) 

22 



<5$ 



(o) 
21 



(o) 
11 



fag) = - V >$<? + (27r(°> + 2^°))$ 
+ (27r(°) + 2/3(°))<l>g ) -4e( )$g ) 



(o) 
12 



(92) 



(93) 



(94) 



Using $jj = 2<pi<j)j, it can be shown that these two equations do not have any new 
information and can be obtained by combining Eqs. (185]) and (IHHj) . 

In the rest of this section, let us us focus on the expansion of the metric up to 
0(r 2 ). First we need the radial derivatives of the spin coefficients as in Eqs. ( 1631) . Let 
us rewrite Eqs. f|63|) keeping this time the contributions of and <p2 & t the horizon, 
and focusing on the spin coefficients which is all we need for expanding the metric up 
to second order: 



,(i) 



(i) 

(0) 







e (1) + 
A« = 



7T 



a 



(i) 
(i) 
(i) 



2 |tt (0) P + 2Re[# 
(o)^(o) _/3(°)) 
(o) 

2 



€ x ~> = 7P '{a 

2//(°)A(°) + vD 4 



f]+4|0i O) | 2 , 
(0)^(0) _ /3 (0) ) 



7P 



(^)) 2 +|A(°)| +2||#| 2 , 
a (i) + ^(i) = 7r (o) /i (o) + 7r (o) A (o ) + v& (o) 

^(D = ^(0)^(0) _^(o) ) + A (o )(/3 (o) 



+ 



(0) 



+ *f + 20( O) 0S O) 



(95a) 
(95b) 
(95c) 
(95d) 
(95e) 
(95f) 
(95g) 
(95h) 



We now need the radial derivatives for the frame functions given in Eqs. 
leads us to the expansion of the frame functions: 



Ol) which 



' U = rk + r 2 (2 |tt(°)| 2 + Re[^ 0) ] + 2|0i O) | 2 ) + 0(r 3 ) , 
n = nfW + r 2 f^TfW + A^ttW + + ^ffi) + 0(r 3 ) 



(0) 



l +r/i (0) +r 2 (/i (O) )2 + |A (O)|2 + | (O) |2 



(0) 



+ 



vY r,l + /- 2 (2^ >A< > + i*£ 0) 



(96a) 

(96b) 
(96c) 

(96d) 
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It is instructive to compare this with Eqs. ( I69p . We see that the electromagnetic field 
modifies ^2 and ^3 according to 

Re[*f ] -> Re[*<°>] + 2|0S O) | 2 , *f -> *<°> + 20i o) 0f . (97) 

In addition, the (/i*- - 1 ) 2 + lA^ 2 term gets an additive correction by \(f>^ | 2 . 

We can finally write the expansion of the metric. It is clear that there is no change 
up to 0(r). Just like the Weyl tensor, the effects of the electromagnetic field first 
enter the metric at second order. Thus, Eqs. ( 1781) and ( !78l) are unchanged and Eq. ( |80l ) 
becomes: 

9% = 4 (k<°>| a - Re[*f ] + 2|<)| 2 ) d {a vd h) v 

+ 4 L(°¥°) + A<«V°> - «f - 2«g o) 0? } ) mg? fyi; (98) 

+4((^) 2 + |A(°)| 2 +i0f i 2 ) m |:» m ;» 

+ (4^°)A(°) - 2*f) m^m® + ^4/i (0) A (0) - 2$ f) fh^fh^ . 

9. Conclusions 

In this paper we have investigated the geometry in the vicinity of a non-extremal weakly 
isolated horizon using the Newman- Penrose formalism. We allow the horizon to have 
generic values of mass, spin and charge and higher multipole moments. We have obtained 
the connection, metric and the Weyl tensor up to 0(r 2 ). On the horizon, we choose 
a parametrization of the null generators so that we get a specific value of the surface 
gravity. In addition, on an initial spherical cross section Sq of the horizon, we specify 
the spin coefficients it, a — A. The spin coefficients 7r and a — (3 out to be time 
independent and they specify respectively the Im[^ 2 ] an d the scalar curvature of S , and 
thus the spin and mass multipole moments of the horizon respectively. The transversal 
expansion and shear fi, A turn out to be time dependent and they determine ^3. In 
addition, we need to specify ^4 on the past-outgoing light cone Af originating from 
Sq. If we consider charged horizons, then we must specify in addition the Maxwell field 
component <\>\ on Sq (this determines the electric and magnetic charges of the horizon) 
and (f>2 on Ao- 

We emphasize once again that this is a local calculation and not all solutions 
produced by this construction are physically meaningful; the interesting issue of global 
existence of solutions has not been addressed here. For example, the uniqueness 
theorems would imply that not all values of the horizon multipole moments would 
lead to asymptotically flat solutions. However, the set of solution we have obtained 
encompasses all physically interesting situations (in so far as the assumption of the 
horizon being isolated is valid). 

There are a number of possible extensions of this work. The first is the study of 
tidally deformed black holes which has, among other things, implications for the the 
equation of motion and self force in general relativity [35]. This requires us to consider 
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Kerr black holes perturbed by a background curvature, and to interpret the resulting 
metric in terms of a moving point particle. This will be studied in a forthcoming paper. 
The near horizon geometry computed here could be matched to an appropriate far zone 
metric to obtain initial data for numerical simulations of binary black hole systems (see 
e.g. [36J) and finally, future gravitational wave observatories might be able to measure 
the tidal deformations through observations of signals from extreme mass ratio systems. 
Knowledge of the near horizon geometry including deviations from the Kerr multipoles 
values (which are allowed here) will be useful in the searches for these signals and in 
testing the Kerr nature of the black hole. Finally beyond the specific calculation of the 
near horizon metric, the characteristic initial value formulation based on the horizon 
might be useful for wave extraction in numerical relativity signals. If So is taken to be a 
cross-section of the horizon at very late times, then the light cone Ao can be viewed as 
an approximate null infinity. The usefulness of this construction remains to be explored. 

Acknowledgments 

I thank Abhay Ashtekar for valuable suggestions and discussions. I am also grateful 
to Bernard Schutz, Toshifumi Futamase, Dirk Piitzfeld and Stas Babak for discussions 
regarding possible applications of this work. 

References 

[1] Ashtekar A, Beetle C, Fairhurst S. Isolated horizons: A generalization of black hole mechanics. 

Class Quant Grav. 1999;16:L1-L7. 
[2] Ashtekar A, Beetle C, Fairhurst S. Mechanics of Isolated Horizons. Class Quant Grav. 

2000;17:253-298. 

[3] Ashtekar A, Fairhurst S, Krishnan B. Isolated horizons: Hamiltonian evolution and the first law. 

Phys Rev. 2000;D62:104025. 
[4] Ashtekar A, et al. Isolated horizons and their applications. Phys Rev Lett. 2000;85:3564-3567. 
[5] Ashtekar A, Beetle C, Lewandowski J. Mechanics of Rotating Isolated Horizons. Phys Rev. 

2001;D64:044016. 

[6] Ashtekar A, Beetle C, Lewandowski J. Geometry of Generic Isolated Horizons. Class Quant Grav. 
2002;19:1195-1225. 

[7] Adamo TM, Newman ET. Vacuum non-expanding horizons and shear-free null geodesic 

congruences. Class Quant Grav. 2009;26:235012. 
[8] Hayward SA. General laws of black hole dynamics. PhysRev. 1994;D49:6467-6474. 
[9] Ashtekar A, Krishnan B. Dynamical horizons: Energy, angular momentum, fluxes and balance 

laws. Phys Rev Lett. 2002;89:261101. 
[10] Ashtekar A, Krishnan B. Isolated and dynamical horizons and their applications. Living Rev Rel. 

2004;7:10. 

[11] Booth I. Black hole boundaries. Can J Phys. 2005;83:1073-1099. 

[12] Gourgoulhon E, Jaramillo JL. A 3+1 perspective on null hypersurfaces and isolated horizons. 

Phys Rept. 2006;423:159-294. 
[13] Hayward SA. Dynamics of black holes; 2008. 

[14] Comeau S, Poisson E. Tidal interaction of a small black hole in the field of a large Kerr black 
hole. PhysRev. 2009;D80:087501. 



The spacetime in the neighborhood of a general isolated black hole 



23 



[15] Poisson E, Vlasov I. Geometry and dynamics of a tidally deformed black hole. Phys Rev. 
2010;D81:024029. 

[16] Poisson E. Metric of a tidally distorted, nonrotating black hole. PhysRevLett. 2005;94:161103. 
[17] Yunes N, Gonzalez JA. Metric of a tidally perturbed spinning black hole. PhysRev. 
2006;D73:024010. 

[18] Damour T, Lecian OM. On the gravitational polarizability of black holes. PhysRev. 
2009;D80:044017. 

[19] Bondi H, van der Burg MGJ, Metzner AWK. Gravitational waves in general relativity. 7. Waves 

from axisymmetric isolated systems. Proc Roy Soc Lond. 1962;A269:21-52. 
[20] Friedrich H. On the Regular and the Asymptotic Characteristic Initial Value Problem for Einstein's 

Vacuum Field Equations. Proc Roy Soc Lond. 1981;A375:169-184. 
[21] Rendall AD. Reduction of the characteristic initial value problem to the Cauchy problem and its 

applications to the Einstein equations. Proc Roy Soc Lond. 1990;427:221-239. 
[22] Friedrich H, Rendall AD. The Cauchy problem for the Einstein equations. LectNotes Phys. 

2000;540:127-224. Einstein's Field Equations and their Physical Interpretation, ed. by B. G. 

Schmidt, Springer, Berlin, 2000. 
[23] Stewart J. Advanced General Relativity. Cambridge Monographs on Mathematical Physics; 1991. 
[24] Lewandowski J. Spacetimes Admitting Isolated Horizons. Class Quant Grav. 2000;17:L53-L59. 
[25] Newman E, Penrose R. An Approach to gravitational radiation by a method of spin coefficients. 

J Math Phys. 1962;3:566-578. 
[26] Diencr P;. Private communication. 

[27] Newman ET, Unti TWJ. Behavior of asymptotically flat empty space. J Math Phys. 1962;3:891- 
901. 

[28] Penrose R, Rindler W. Spinors and Stacetime: 1. Two-Spinor Calculus and Relativistic Fields. 
Cambridge, Uk: Univ. Pr. ( 1984) 458 P. ( Cambridge Monographs On Mathematical Physics); 
1984. 

[29] Chandrasckhar S. The mathematical theory of black holes. Oxford Classic Texts in the Physical 
Sciences; 1985. 

[30] Goldberg JN, MacFarlane AJ, Newman ET, Rohrlich F, Sudarshan ECG. Spin s spherical 

harmonics and cdth. J Math Phys. 1967;8:2155. 
[31] Dreyer O, Krishnan B, Shoemaker D, Schnetter E. Introduction to Isolated Horizons in Numerical 

Relativity. Phys Rev. 2003;D67:024018. 
[32] Beetle C. Approximate Killing Fields as an Eigenvalue Problem; 2008. 

[33] Cook GB, Whiting BF. Approximate Killing Vectors on S**2. PhysRev. 2007;D76:041501. 

[34] Ashtekar A, Engle J, Pawlowski T, Van Den Broeck C. Multipole moments of isolated horizons. 

Class Quant Grav. 2004;21:2549-2570. 
[35] Poisson E. The motion of point particles in curved spacetime. Living Rev Rel. 2004;7:6. 
[36] Yunes N, Tichy W. Improved initial data for black hole binaries by asymptotic matching of 

post-Newtonian and perturbed black hole solutions. PhysRev. 2006;D74:064013. 



